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Abstract 

It is shown that under a certain condition on a semimartingale and 
a time-change, any stochastic integral driven by the time-changed semi- 
martingale is a time-changed stochastic integral driven by the original 
semimartingale. As a direct consequence, a specialized form of the Ito 
formula is derived. When a standard Brownian motion is the original 
semimartingale, classical Ito stochastic differential equations driven by 
the Brownian motion with drift extend to a larger class of stochastic dif- 
ferential equations involving a time-change with continuous paths. A form 
of the general solution of linear equations in this new class is established, 
followed by consideration of some examples analogous to the classical 
equations. Through these examples, each coefficient of the stochastic dif- 
ferential equations in the new class is given meaning. The new feature 
is the coexistence of a usual drift term along with a term related to the 
time-change. 



1 Introduction 

Among the most important results in the theory of stochastic integration is 
the celebrated Ito formula, which establishes a stochastic calculus for stochastic 
integrals driven by a semimartingale. In general, given a d-dimcnsional semi- 
martingale X = (X 1 , . . .,X d ) starting at 0, if / : R d — > R is a C 2 function, 
then f(X) is a one-dimensional semimartingale, and, for all t > 0, with proba- 
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One useful implication of the Ito formula (1.1) is the product rule. Namely, if 
Y and Z are both one-dimensional scmimartingales starting at 0, then, for all 
t > 0, with probability one 



These formulas arc indispensable tools for working with stochastic differential 
equations (SDEs). 

Our motivation to investigate stochastic integrals driven by a time-changed 
semimartingalc originated in a desire to develop a stochastic calculus when the 
time-change is the first hitting time process of a stable subordinator of index 
between and 1. Meerschaert and Schcffler [10, 11] show that this type of 
process arises as the scaling limit of continuous time random walks. 

Section 2 first introduces the significant concept of synchronization, which 
connects a semimartingale with a time-change in a certain manner. A time- 
change (T t ) is a cadlag, nondecrcasing family of stopping times. Given a one- 
dimensional semimartingale Z starting at 0, the composition of Z and T, de- 
noted ZoT or (ZT t ), is called the time-changed semimartingale. We occasionally 
refer to t and T t as the original clock and the new clock, respectively. With the 
notion of synchronization. Jacod [6] explains how to recognize a time-changed 
stochastic integral of the form J Q Tt H s dZ s in terms of an integral with respect to 
the time-changed semimartingale (Zr t ) (Lemma 2.3). However, this statement 
does not answer the following question: 

Q: When and how can a stochastic integral J* K s dZT 3 driven by a time- 
changed semimartingale be realized by way of an integral driven by the 
original semimartingale (Z t )7 

In Section 3, Theorem 3.1 provides a complete answer to the above question. 
Namely, K s dZT 3 = Jq* Kgi s _\dZ s , where S is the first hitting time process of 
T. An important corollary of Theorem 3.1 is a form of the Ito formula (1.1) for 
a C 2 function of a process which contains a stochastic integral driven by a time- 
changed semimartingale (-ZsJ where (E t ) is a continuous time-change, meaning 
a time-change with continuous paths. The formula can be reexpressed in terms 
of usual stochastic integrals driven by the original semimartingale and the con- 
tinuous part of the semimartingale's quadratic variation. A generalization of 
this formula is a time-changed Ito formula provided in Theorem 3.3. 
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Theorem 3.1, from which the time-changed Ito formula is derived, can be 
regarded as a powerful tool in handling a new class of SDEs which are driven by 
Lebesgue measure, a continuous time-change, and a time-changed scmimartin- 
gale (Section 4). The simplest, yet quite significant subclass, of such SDEs are 
ones with linear coefficients: 

(1.3) dX t = (p 1 (t,E t ) + p 2 (t,E t )X t )dt+(p 1 (t,E t ) + p 2 (t,E t )X t )dE t 

+ (<j 1 (jk,Et) + a 2 (t,Et)Xt)dBi3 t , 

where B is a standard Brownian motion. The new feature of this class of SDEs is 
the coexistence of a term representing a drift under the new clock E t along with 
a usual drift based on the original clock t. Theorem 4.5 establishes a general 
form of the solution to SDE (1.3), in which again Theorem 3.1 is applied to 
obtain another representation of the solution. 

Section 5 compares some SDEs of the form (1.3) with classical Ito SDEs, 
described as 

(1.4) dY t = (h(t) + b 2 (t)Y t )dt + (r^t) + T 2 (t)Y t )dB t . 

The comparison reveals the role of the dE t term appearing in SDE (1.3). Namely, 
fij can be ascribed to either bj or tj in (1.4), depending on the way the model 
(1.3) is constructed (Remark 4.7 (b)). These examples also illustrate methods 
for obtaining statistical data of the solution, such as the mean and variance. 

2 Preliminaries — Stochastic Integrals and 
Time-changes 

Throughout this paper, a complete filtered probability space (f2, !F, {Ft), P) is 
fixed, where the filtration (J- t ) satisfies the usual conditions; that is, it is right- 
continuous and contains all the P-null sets in T . For simplicity, unless mentioned 
otherwise, all processes are assumed to take values in R and start at 0. 

A process Z is said to be cadlag (resp. cdgldd) if Z has right-continuous 
sample paths with left limits (resp. left-continuous sample paths with right 
limits). The assumption that Z is cadlag or caglad requires the process to have 
at most countably many jumps. Associated to a cadlag process Z is its jump 
process AZ = (AZt) where AZt '■= Zt — Zt- with Zt_ denoting the left limit at 
t and Zo- = by convention. Let B(J-t) and h(!Ft) respectively denote the class 
of cadlag, (J-t)-adapted processes and that of caglad, (J-t)-adapted processes. 

A cadlag process Z is an (Tt)-semimartingale if there exist an (^ r t )-local 
martingale M and an (J r t )-adapted process A of finite variation on compact sets 
such that Z = M+A. Although this decomposition is not unique in general, the 
local martingale part M can be uniquely decomposed as M = M c + M d with 
a continuous local martingale M c and a purely discontinuous local martingale 
M d . The process M c is determined independently of the initial decomposition of 
Z into M and A, and Z c is defined to be the unique continuous local martingale 
part M c of Z; i.e., Z c := M c ([7, I. Prop. 4.27]). 
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The class of semimartingales forms a real vector space which is closed un- 
der multiplication. It is known to be the largest class of processes for which 
the Ito-type stochastic integrals arc defined. The notion of predictability is es- 
sential in the construction of a stochastic integral driven by a semimartingalc. 
Let V(Tt) be the smallest c-algebra on R + x f2 which makes all processes in 
L(^ r t ) measurable. An (Tt)- predictable process is a process which is V(Tt)- 
measurablc. Given an (^-" t )-semimartingalc Z, let L(Z,Tt) denote the class of 
(^)-predictable processes H for which a stochastic integral driven by Z, denoted 
(H • Z) t = L H s dZ s , can be constructed. A brief summary of the construction 
appears in Appendix. 

The quadratic variation of a semimartingalc Z, denoted [Z,Z], can be de- 
fined by way of a stochastic integral. It is the cadlag, (J r t )-adapted, nondecreas- 
ing process given by 



By polarization, the map [•, •] becomes a symmetric, bilinear form on the class 
of semimartingales. For semimartingales Y and Z, note that [Y, Z] c does not 
denote its continuous martingale part, which is of course zero, but it is defined 
to be its continuous part; namely, 



It follows by comparing this definition with Theorem 4.52 in [7, Chap. I] that 



The following are some of the basic but key properties of stochastic integrals 
which will be employed in the subsequent sections. 

Properties 2.1. Let Y and Z be (Tt)- semimartingales. Let H G L(Z,Tt). 

(1) H • Z is again an (Ft) -semimartingalc 

(2) A(H»Z) = H ■ AZ. 

(3) Additionally, if H G L(Y, T t ), then H • (Z + Y) = H • Z + H »Y . 

(4) IfJe L(H • Z, T t ), then J • (H • Z) = (J ■ H) • Z . 

(5) IfKe L(Y, T t ), then [H • Z, K • Y] = (H ■ K) • [Z, Y] . 

An (Tt)-time- change is a cadlag, nondecrcasing family of (J-" t )-stopping times. 
It is said to be finite if each stopping time is finite almost surely. Let (T t ) be a 
finite (.F^-time-change and define a new filtration (Q t ) by Q t = Tt ± - Then (Q t ) 
also satisfies the usual conditions since the right-continuity of (Tt) and (T t ) im- 
plies that of (Qt)- In addition, for any (F" t )-adapted process Z, the time-changed 
process (Zr t ) is known to be (C/ t )-adapted. In fact, more can be said. 



(2.1) 




[Y, Z] c t := [Y, Z] t - ]T A[Y, Z) s = [Y, Z] t — ^ AY S ■ AZ S . 




[Y,Z] C = [Y C ,Z C ]. 
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Lemma 2.2. ([6, Cor. 10.12]) Let Z be an (Tt)- semimartingale. Let (T t ) be a 
finite (Ft) -time- change. Then (ZT t ) is a (Q t )- semimartingale where Q t '■= Tt^ 

Furthermore, if Z is an (.F t )-adapted process of finite variation on compact 
sets, then (Z-rJ also has finite variation on compact sets. On the other hand, 
the local martingale property, in general, is not preserved under a given time- 
change. In other words, even if Z is an (J r t )-local martingale, the time-changed 
process (Zx t ) may fail to be a (C/ t )-local martingale. A simple example is a 
standard ( Tt )-Brownian motion Z = B with the finite (J-t)-time-change (T t ) 
defined by T t := inf{s > 0; B s = t}. In this case, Bt ± = t for every t > 0. Thus, 
the time-changed Brownian motion is no longer a local martingale. 

One way to exclude this unexpected possibility is to introduce the notion 
of synchronization, which turns out to be an essential concept in developing a 
stochastic calculus for integrals driven by a time-changed semimartingale. A 
process Z is said to be in synchronization with the time-change (T t ) if Z is con- 
stant on every interval [T t _, T t ] almost surely. We occasionally write Z ^ S ynch T 
for shorthand. Other properties that a time-change preserves appear in [6, Thm. 
10.16]. In the literature, Jacod [6], Kallsen and Shiryaev [8] use the expression 
(Tt)- adapted in describing a process being in synchronization with a time-change 
(Tt). A different terminology (Tt)- continuous is used by Revuz and Yor [14]. 
Nevertheless, the term synchronization is adopted here to avoid any possible 
confusions or misunderstandings that the other expressions may create. 

One quite simple yet significant observation, which connects the notion of 
synchronization with stochastic integrals, is that if an (.F t )-semimartingale Z 
is in synchronization with a finite (.F t )-time-change (T t ) and if H G L(Z,T t ), 
then (H T ( t -)) <= L(Z o T,Q t ), where, H T u_\ denotes the process H evaluated 
at the left limit point T t - of T at t. This observation leads to the consid- 
eration of two integral processes (J Q H s dZ s ) and (J Q H T ^ s _)dZT s ). By Prop- 
erty 2.1 (1), these are semimartingalcs with respect to the nitrations (Tt) and 
(Qt), respectively. By Lemma 2.2, the former stochastic integral can be time- 
changed by (T t ) to produce another (C^)-semimartingalc. Jacod [6] shows that 
the two (C/ t )-semimartingales (J^* H s dZ s ) and (J* H T ^ s _^dZT 3 ) coincide for any 
H G L(Z,Tt). This fact plays a significant role in establishing the basic Theo- 
rem 3.1; hence, it is stated here as a lemma. 

Lemma 2.3. (1st Change-of- Variable Formula [6, Prop. 10.21]) Let Z 

be an (Tt) -semimartingale which is in synchronization with a finite (Tt)-time- 
change (T t ). Lf H G L(Z,Tt), then (H T ( t _)) <= L(Z oT,Q t ) where Qt '■= Tr t - 
Moreover, with probability one, for all t > 0, 



Lemma 2.4. ([6, Thm. 10.17]) Let Z be an (Tt)- semimartingale which is 
in synchronization with a finite (Tt) -time- change (T t ). Then Z c and \Z,Z\ are 
also in synchronization with (T t ). Moreover, 



(2.2) 




(2.3) 



[ZoT,ZoT] = [Z,Z]oT, (ZoT) c = Z c oT. 
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The following simple example explains the significance of the synchronization 
assumption in Lemmas 2.3 and 2.4. 

Example 2.5. Let Z = B be a standard (jF t )-Brownian motion, and define a 
deterministic time-change (T t ) by T t ■= Iri iOC )(t), where Ia denotes the in- 
dicator function over a set A. Let H be a deterministic process given by 
Ht = I(i/2, cx>)(*), then H T{t _ } = I (lj00 )(t). Hence, 

Ti rl rl 

H s dB s = / H s dB s = / dB s =B 1 -B 1/2 ; 

JO Jl/ 2 

H T{s „ } dB Ts = [ (iB Ts = 0. 
Jo 

Therefore, the two integrals in (2.2) fail to coincide. Moreover, it follows from 
(2.1) that 

[B o T, B o T]i = (73 Tl ) 2 - 2 / B Ts _dB Ts = B\-2 ( dB T> = S^, 

Jo Jo 

whereas the fact that [B,B] t = t yields {[B,B] o T)i = Ti = 1. Therefore, the 
first equality in (2.3) does not hold. Furthermore, since BoT is not a continuous 
process, (BoT) c and B c oT(= BoT) fail to coincide. Thus, the second equality 
in (2.3) does not hold either. Note that the Brownian motion B never stays 
flat on any time interval, and hence is not in synchronization with the above 
time- change (Tt). □ 

The next lemma will be used in the proof of Theorem 3.1. 

Lemma 2.6. Let Z be an (Tt)-semimartingale which is in synchronization with 
a finite (!Ft) -time-change (Tt). Let H € L(Z,J- t ). Then the stochastic integral 
hi • Z is also in synchronization with (Tt). 

Proof. Fix t > 0, and let u£ [Tt-, T t ] . Since Z ^ syn ch T, Z is constant on [it, T t ] ; 
hence, (H • Z) Tt - (H • Z) u = J^H s dZ s = 0. Therefore, (H»Z) Tt = (H»Z) U . 
Thus, H • Z is constant on [T t _, T t ] . □ 

The following lemma and its corollary clarify the situation of main concern 
in this paper. The first hitting time process of a given cadlag, nondecreasing 
process S is a process T defined by T t := inf{it > 0; S u > t}. It is easy to see 
that T is also cadlag and nondecreasing. Note that every (y r t )-adapted, cadlag, 
nondecreasing process has paths of finite variation on compact sets; hence, a 
priori it is an (y r t )-semimartingale. 

Lemma 2.7. 

(1) Let S be a nondecreasing (Tt)-semimartingale such that Imit—^oo St = oo 
and St = only when t = 0. Then the first hitting time process T of S is a 
finite (Tt)-time change such that lim^oo T t = oo. Moreover, if S is strictly 
increasing, then T has continuous paths. 
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(2) Let T be a finite {Ft) -time- change such that lim^oo T t = oo and T t — 
only when t = 0. Then the first hitting time process S of T is a nonde- 
creasing [Tt)-semimartingale such that lim^oo St = oo. Moreover, if T 
has continuous paths, then S is strictly increasing. 

Proof. (1) The assumption that St > for t > and lim^oo St = oo implies 
that To = and each random variable T t is finite. In addition, since each St is 
a real-valued random variable, it follows lim^oo T t = oo. Fix t > 0. Since S is 
(.F t )-adapted, {T t < s} = {S s - > t} G .Ts- C .F s for any s > 0, and obviously 
{T t < 0} = G JF . Hence, T t is an (jF t )-optional time. It follows from the 
right-continuity of (T t ) that T t is an (.T^-stopping time. (See [9, Prop. 1.2.3].) 
Therefore, T is a finite (jT t )-time change. Moreover, if S is strictly increasing, 
then T obviously has continuous paths. 

(2) The assumption that T t > for t > and lim^oo T t = oo implies 
that Sq = and each random variable St is finite. In addition, since each T t 
is a real-valued random variable, it follows lim^oo = oo. Fix s > 0. For 
any t > 0, since T t _ is also an (jF t )-stopping time, {S s > t] = {T f _ < s} G 
T s . Also, {5 S > 0} = il G .F s . Hence, 5 S is JFg-measurable. Therefore, 
S is ( )-adapted. Since S is also cadlag and nondecreasing, it is an (J-*)- 
semimartingale. Moreover, if T has continuous paths, then it is clear that S is 
strictly increasing. □ 

Remarks 2.8. 

(a) Lemma 2.7 establishes that a nondecreasing (J-t)-semimartingale S and a 
finite (^-"t)-time-change T are 'dual' in the sense that cither process with the 
specified condition induces the other. 

(b) Part (1) of Lemma 2.7 assumes that lim^oo St = oo and St = only when 
t = 0. The first condition ensures that T does not blow up in finite time. We 
may lift this condition by restricting attention to Tt with t G [0,i*) where 
£, = sup 0<s<oo S s , the explosion time of T. However, the condition St = 
only when t = cannot be weakened. It guarantees that the first hitting 
time process T starts at 0. The same argument applies to the assumption 
on T in Part (2). 

Notation 2.9. In light of Remark 2.8 (a), for a pair of a nondecreasing 
(^-t)-semimartingale S and a finite (Tt )-time-change T, [S i — ► T] and 
[S < — i T] are used to indicate respectively that 5* induces T and that T induces 
S as described in Lemma 2.7. If S is strictly increasing and T has continuous 
paths, then the double brackets [[S i — > T]] and [[5 < — i T}] are employed instead. 
Hence, the double bracket notation assumes stronger conditions than the single 
bracket notation. Hereafter, the notation D and E will be used to denote a 
pair of a strictly increasing semimartingale and a continuous time-change. This 
notation is chosen to be compatible with the continuous time-change E, which 
is induced by a strictly increasing, stable subordinator D of index between 
and 1, in the papers of Meerschaert and Scheffler [10, 11] on continuous time 
random walks. 
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3 Stochastic Calculus for Integrals Driven by a 
Time-changed Semimartingale 

This section establishes a stochastic calculus for integrals driven by a time- 
changed semimartingale. The central problem is to understand such integrals 
by rephrasing them in terms of integrals driven by the original semimartingale. 
Solving this problem is almost equivalent to providing a way to recognize SDEs 
driven by a time-changed semimartingale, which aids the analysis of problems 
that appear in applications. 

The following theorem, at first glance, may seem quite simple, but its im- 
pact on the formulation of our stochastic calculus is profound. Recall that all 
processes, unless specified otherwise, are assumed to take values in K and start 
at throughout the paper. 

Theorem 3.1. (2nd Change-of- Variable Formula) Let Z be an (J 7 t)- 
semimartingale. Let S and T be a pair satisfying [S i — ► T] or [S < — i T] . 
Suppose Z is in synchronization with T. Lf K G L(Z o T,Qt), then (-Kg( t _)) G 
L(Z,Qs t ) where Gt : = ^T t - Moreover, with probability one, for all t>0, 

(3.1) / K s dZ Ts = / * K S (s-)dZ s . 
Jo Jo 

Proof. By Lemma 2.2, both T and Y '■= ZoT arc (tj t )-scmimartingales. Since T 
is a nondecreasing (^t)-scmimartingale such that limt—xxj Tt = oo and To = 0, it 
follows from Part (1) of Lemma 2.7 along with Remark 2.8 (b) that S is a finite 
(St)-time-changc. On any half open interval [S S -,S S ), T is obviously constant 
by construction and hence so is Y . Moreover, since Z ~ sync h T, 

(Z o T) S(s ) = Z T(S ( S )) = Z T(s(s) _) = Z T (s( a -)) = (Z o T) S ( S _). 

Hence, Ys s = Ygu_\. Thus, Y is constant on any closed interval [5 S _,5 S ]. 
Therefore, Y ~ sync h S. 

Now, let K G L(Y,Q t ). Then it follows from Lemma 2.3 that (K S ( t _)) G 
L(Y o S,Gs t )- By the 1st change-of-variable formula (2.2) and the assumption 
Z ~ S ynch T, with probability one 

St pt pt pt 

K s dY s = I K s{s ^ ) dY Ss = l K s ^_)dZ T ^ s)) = I K S ( s -)dZ s 
'o Jo Jo Jo 

for alH > 0. Hence, with probability one, 

(3.2) [ STt K s dY s = Ks {a -)dZ, 
Jo Jo 

for all t > 0. Since Y ^ S ynch S, Lemma 2.6 yields K • Y ~ S ynch S. Any t is 
contained in the interval [5r(t)_, 5r t ], so {K • Y)s T — {K •Y) t . Thus, (3.2) 
establishes (3.1). □ 



Remarks 3.2. 



(a) Theorem 3.1 guarantees that any stochastic integral driven by a time- 
changed semimartingale is a time-changed stochastic integral driven by the 
original semimartingale, as long as the semimartingale is in synchronization 
with the time-change. 

(b) If a pair D and E satisfies \[D \ — > E]) or [[D < — i E]], then any process 
Z is automatically in synchronization with E due to the continuity of E. 
Therefore, under cither of these stronger conditions, Theorem 3.1 is valid 
for an arbitrary (J r t )-semimartingale Z . 

In light of Remark 3.2 (b), when [[£) i — ► E}} or [[D < — i E]}, the Ito formula 
for stochastic integrals driven by a time-changed semimartingale can be refor- 
mulated in a nice way via the 2nd change-of- variable formula (3.1) obtained in 
Theorem 3.1. The proof of Theorem 3.3 is provided in full detail since it demon- 
strates important computational techniques on quadratic variations which are 
frequently employed in Section 4. 

Theorem 3.3. (Time-changed Ito Formula) Let Z be an (Ft) -semimartingale. 
Let D and E be a pair satisfying [[D i — > E\] or [[D < — i E\\. Define a filtration 
{Qt) by Q t := Te ± ■ Let X be a process defined by 

(3.3) X t : = (A • m) t + (F • E) t + (G • {Z o E)) t 

A s ds+ [ F s dE s + [ G s dZ Es 



where A € L(m, Q t ), F € L(E, Q t ), G G L(Z o E, Q t ), and m is the identity map 
on K corresponding to Lebesgue measure. If f : M — > R is a C 2 function, then 
f(X) is a (Qt) -semimartingale, and with probability one, for all t > 0, 

(3.4) f(X t )-f(0)= [ f'(X s _)A s ds+ [ '/'(^(a-)-)^^-)^ 
Jo Jo 

" f'{X D{s ^)G D{s ^dZ s + X - f"(X D(s ^){G D{s ^} 2 d[Z, Z] c s 



+ {f(Xs) - f(X s -) - f'(X s -)AX s }. 

0<s<t 

In particular, if Z is a standard Brownian motion B, then with probability one, 
for all t > 0, 

(3.5) f(X t )-f(0)= / f'(X s )A s ds+ / f'{X D{s ^)F D{s ^d S 

JO JQ 

+ r/'(^ (s -))G D(H rfB s + J ! ' f"{X D{s _ ) ){G D{s _ ) } 2 ds. 
Jo 1 Jo 
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Proof. Since the process X in (3.3) is defined to be a sum of stochastic inte- 
grals driven by (C/ t )-semimartingales, X itself is also a (C/t)-semimartingale by 
Property 2.1 (1). The Ito formula (1.1) with d = 1 yields, for all t > 0, 

(3.6) f(X t )-f(0)= f f'{X^)dX s + \[ f"(X s .)d[X,X] c a 

Jo 1 Jo 

+ £ {f(X s )-f(X s .)-f'(X s -)AX s }. 

0<s<t 

Using Properties 2.1 (3), (4), 

(3.7) / f(X s _)dX s =[ f'(X s _)A s ds+ f f'(X s _)F s dE s 
Jo Jo Jo 

+ f f'(X s _)G s dZ Es . 
Jo 

By Remark 3.2 (b), any semimartingale is in synchronization with E. Hence, 
the 2nd changc-of-variable formula (3.1) yields 

(3.8) f f'(X a -)dX a =f f'(X a -)A a ds+ [ '/'(V)-)^!-)^ 
Jo Jo Jo 

/•Et 

+ / f (XD{s-)-)GD{s-)dZ s . 

Jo 

For the second integral on the right hand side of (3.6), first let Y := Z o E. 
We claim that 

(3.9) [X,X]° t = [ G 2 s d[Y,Y] c s . 

Jo 

To prove this, first note that m and E are both continuous processes of finite 
variation on compact sets. By [13, II. Thm. 26], 

[m,Y] t = Y, A K y ]«= (Am s ) ■ (AX S ) = 

0<s<t 0<s<t 

for all t > 0. Hence, [to, Y] = 0. Similarly, [m, m] = [m, E] = [E, E] = [E, Y] 
0. Therefore, the bilinearity of [•, •] and Property 2.1 (5) imply 

(3.10) [X, X] = [A • to + F • E + G • Y , A»m + F»E + G»Y] 

= G 2 .[y,y]. 

Now, let J t := J2o<s<t Y U so that [Y,Y] c t = [Y, Y] t ~ J t . Then the pure 
jump process, J, shares with [Y, Y] the same jump times and sizes. Therefore, 

GlA[Y,Y] s = G 2 AJ S = fcldJ s . 

0<s<t 0<s<t J ° 
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Hence, it follows from (3.10) together with Properties 2.1 (2), (3) that 
[X,X]° = [X,X} t - £ A[X,X] s = (G 2 .[Y,Y)) t - £ G 2 S A[Y, Y] s 

0<s<t 0<s<t 

= f G 2 3 d[Y,Y] s - f G 2 s dJ s = f G 2 s d[Y,Y]% 
Jo Jo Jo 

thereby establishing (3.9). 

Since Z ~ S ynch E, repeated use of Lemma 2.4 yields 

(3.11) [Y, Y] c = [Y c , Y c ] = [Z c oE,Z c oE} = [Z c , Z c ]oE= [Z, Z] c o E. 
Together (3.9) and (3.11) yield 

[X,X}<= f Gld[Z,Z] c Bs . 
Jo 

Therefore, it follows from Property 2.1 (4) and the 2nd change-of- variable for- 
mula (3.1) that 

(3.12) f f"(X s .)d[X, Xf 8 = f f"(X s „)G 2 s d[Z, Z\l s 
Jo Jo 

Et r{x D{s ^){G D{s ^} 2 d[z,zr s . 

Equality (3.4) follows by plugging (3.8) and (3.12) into Formula (3.6). 

If Z = B is a standard Brownian motion, then the continuity of m, E 
and B o E together with Property 2.1 (2) imply X is also continuous. Since 
[B,B\1 = [B,B] t = t, statement (3.5) follows immediately. □ 

A similar proof yields the multidimensional version of Theorem 3.3. For a 
multidimensional process W, its i-th component is denoted W z . 

Corollary 3.4. Let Z be an n- dimensional {Tt)-semimartingale starting at 0. 
Let D and E be a pair satisfying [[D i — > E]] or [[D < — i E]]. Define a filtration 
(Gt) by Qt '■= J~E t - Let X be a d- dimensional process defined by 

X t := / A s ds + / F s dE s + V / G k s dZ k Es 
Jo Jo k=l 

where A, F and G k — (G k ' 1 , . . . , G k ' d ) (k ~ 1, . . . ,n) are d- dimensional pro- 
cesses for which all the above integrals are defined. If f : Mr — > M is a C 2 
function, then f(X) is a {Q t )-semimartingale, and with probability one, for all 
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t > 0, 

(3.13) f(X t )-f(0) 



d " f E t as 

+EE / ^(^(.-)-)^Jo^ 

i,j = l k,£=l JV 



a/ 

0<s<t " i=l 

Remarks 3.5. 

(a) The first integral in Formula (3.4) can also be expressed as a time-changed 
stochastic integral. By the 2nd change-of- variable formula (3.1), 

(3.14) f f'{X s _)A s ds 
Jo 

= f f(X s _)A 8 dD Es + Y, f(X 3 -)A s A(DoE) s 

- f Et f{X D(s _ ) _)A D(s _ ) dD s + V f(X s _)A s A(DoE) s 

Jo 0<s<t 

as long as all integrals are defined. The additional term arises due to the 
discontinuities of D. 

(b) The stronger condition [[D i — > E]] or \[D < — i E}], rather than [D i — > 
E] or [D < — i E], is essential in establishing the nice representations (3.4) 
and (3.13). For example, if E has jumps, then the stochastic integral 

f'(X s _)F s dE s in (3.7) may not be rephrased as a time-changed integral 
driven by ds since the identity map m(s) = s is no longer in synchronization 
with E. Moreover, the equalities [E, E] = and [E, Y] = both may fail, 
which implies more terms need to be included in (3.9). 

(c) In real situations, the distributions of Z, D and E are known through 
statistical data, and scientists will seek to reveal the behavior of a process 
X described via an SDE of the form 

(3.15) dX t =p(t,Et,Xt)dt + n(t,Et,Xt)dEt+a(t,Et,Xt)dZ E f 

Formula (3.4) encourages handling the solution to Equation (3.15) via con- 
ditioning. In particular, when Z is continuous and A = 0, the right hand 
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side of Formula (3.4), conditioned on E t , can be regarded as usual stochastic 
integrals driven simply by Lebesgue measure, Z and [Z,Z]. 

The following example provides a sense of the kinds of results that can be 
obtained using Theorem 3.3 together with conditioning. 

Example 3.6. Let D be an (.F^-stable subordinator of index (3 € (0, 1) which 
is independent of a standard (J-" t )-Brownian motion Z — B. The process D 
is strictly increasing. Let E be the associated continuous time-change so that 
[\D i — > E]]. Then under a certain condition, the transition probability density 
p x (t, x, y) of a solution X to the SDE 

(3.16) dX t = fx(X t )dE t + a(X t )dB Et with X Q = x 

satisfies the following time-fractional partial differential equation (PDE) in the 
weak sense: 

(3.17) n?p x (t, x, y) = - A { f i(y)p x (t , x, y) } + ^ {a 2 (y)p x (t,x,y)}, 

with initial condition p x (0,x,y) = S x (y). Here, D* is the Caputo fractional 
derivative of order (3 with respect to the time variable t (see [4]), and S x is 
the Dirac delta function with mass at x. For the proof, see Hahn, Kobayashi 
and Umarov [5, Thm. 4.1]. Furthermore, that paper provides a more general 
perspective on this matter in the framework of time-changed Levy processes 
and their associated pseudo-differential equations. Moreover, the above result 
is derived there without the use of Theorem 3.3, but based on Theorem 4.2 
of the present paper. The advantage of the approach which employs the time- 
changed ltd formula (3.5) is that it reveals the connection between the stochastic 
calculus for a time-changed Brownian motion and its associated time-fractional 
PDE (3.17). A further remark on Equations (3.16) and (3.17) is provided in 
this paper in Example 5.4. □ 



4 SDEs Including Terms Driven by a Time-changed 
Semimartingale 

A classical Ito SDE is of the form 

dY t = b(t,Y t )dt + T(t,Y t )dB t 

where B is a standard Brownian motion. As stated in Remark 3.5 (c), the 2nd 
change-of- variable formula (3.1) is a useful tool in handling a larger class of 
SDEs of the form 

dX t = p(t, E t ,X t )dt + fifa Et,X t )dEt + a(t, E t ,X t )dB Et , 

where E is a continuous time-change. Note that the sample path t i— ► Et is 
not necessarily absolutely continuous with respect to Lebesgue measure; hence, 



13 



the dE t term appearing above in general cannot be rewritten in terms of dt. 
For example, if E is the first hitting time process of a stable subordinator D of 
index between and 1, then the sample path t 1— > E t is flat almost everywhere. 
Therefore, if E t had a representation E t = J g(s)ds for some integrable function 
g, then it would follow that Et = for all i > 0, contradicting the fact that 
Iirut-too E t = oo. More generally, if E is the first hitting time process of a 
strictly increasing Levy process with infinite jumps and no drift, then E is not 
absolutely continuous with respect to Lebesgue measure. For definition and 
properties of Levy processes, consult [1] or [15]. 

The new feature of this larger class of SDEs is the coexistence of a usual drift 
term along with a term representing a factor ascribed to the time-change. The 
aim of this section is to provide ways of recognizing this new larger class of SDEs 
by analyzing their solutions and making comparisons between the two classes of 
SDEs. For a general treatment of classical Ito SDEs, see [9] or [12]. Regarding 
methods for obtaining explicit forms of solutions to classical Ito SDEs, consult 
[3, Chap. 4]. Many basic models are introduced in [16] with an abundance of 
interpretations and insights. 

Let Z be an (J r t )-scmimartingalc and let E be a continuous (J r t )-time-change. 
The general form of SDEs discussed here is 

(4.1) dX t = p(t,E t ,Xt-)dt + v(t,E u Xt-)dE t + <j(t,E t ,X t -)dZ Ei 

with Xq = xq, 

which is understood in the following integral form: 

(4.2) X t = x Q + [ p{s,E a ,X a -)ds+ f fx(s,E s ,X s _)dE s 

Jo Jo 

+ f o{s,E s ,X s ^)dZ E „ 
Jo 

where xq is a real constant, and p, p, a are real-valued functions, defined on 
K_i_ x M + x K, which satisfy the following Lipschitz condition: there exists a 
positive constant L such that 

(4.3) \p(t,u,x) - p(t,u,y)\ + \fi(t,u,x) - p(t,u,y)\ 

+ \a(t,u,x) - a(t,u,y)\ < L\x - y\ 

for all t,u € M+ and x, y £ M. For technical reasons, we require assumption 

(4.4) X G D(ft) 

(p(t,Et,X t -)), (n{t,Et,X t _)), (*(t,Et,Xt-))eUg t ), 

where Q t := Te ± - One example of such functions is a 'linear' map p(t,u,x) = 
pi(t, u) + pi(t, u) ■ x, where p\, p2 are bounded continuous functions on R + xR + . 

Lemma 4.1. (Existence and Uniqueness of Solution) Let Z be an {Tt)- 
semimartingale. Let D and E be a pair satisfying [[D < — > E]] or [[D < — i E\\. 
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Suppose p, (i, a are real-valued functions defined on K + x R + x M satisfying 
Lipschitz condition (4.3) and assumption (4.4). Then there exists a unique 
{Qt)-semimartingale X for which (4.1) holds, where Gt : = ^e,- 

Proof. The identity map m corresponding to Lebesguc measure can be regarded 
as a (<7 t )-semimartingale. Moreover, E and Z oE are also (£t)-semimartingalcs 
due to Lemma 2.2. The existence and uniqueness of a strong solution X to SDE 
(4.1) is guaranteed by conditions (4.3) and (4.4), upon reformulating Theorem 
7 of [13, Chap. V] with operators Fj : ©(&) — ► L(&) (j = 1, 2, 3) defined by 

F 1 (A-) t = p(t,^ J Jf t _), F 2 (X) t = n(t 1 E t ,X t -), F 3 (X) t = a(t,E t ,X t _). 

Furthermore, it follows from Property 2.1 (1) and the integral expression (4.2) 
that X is a (^t)-scmimartmgale. □ 

Now that the existence and uniqueness of a solution to an SDE of the form 
(4.1) is established, the following two SDEs both make sense: 

(4.5) dX t = (i(Et,X t -)dE t +(T(Et,Xt-)dZE t with X = x ; 

(4.6) dY t = fi(t, Y t -)dt + a(t, Y t J)dZ t with Y = x . 

Together the change-of- variable formulas (2.2) and (3.1) yield Theorem 4.2, 
which in turn reveals a close connection between the classical Ito-type SDE 

(4.6) and our new class of SDEs in (4.5). 

Theorem 4.2. (Duality of SDEs) Let Z be an [Tt)-semimartingale. Let D 
and E be a pair satisfying [[D i — > E\] or [[D < — i E]]. 

(1) If a process Y satisfies SDE (4.6), then X := Y o E satisfies SDE (4.5). 

(2) If a process X satisfies SDE (4.5), then Y := X o D satisfies SDE (4.6). 

Proof. (1) Suppose Y satisfies SDE (4.6), and let X := YoE. Since any process 
is in synchronization with the continuous (J r t )-time-change E, the 1st change- 
of- variable formula (2.2) yields 

(4.7) X t =x + [ * fx{s,Y s _)ds + [ * a{s,Y s _)dZ s 

Jo Jo 




In general, the equality Yb(s)- = {Y o E) s - may fail. The failure can occur 
only when E is constant on some interval [s — e, s] with e > 0. However, the 
integrators E and Z o E on the right hand side of (4.7) are constant on this 
interval; hence, the difference between the two values Y E ^)_ and (YoE) s - docs 
not affect the value of the integrals. Thus, (4.7) can be reexpressed as 

(4.8) X t =x + f fi(E s ,(YoE) s _)dE s + f a{E s , (Y o dZ E „ 
Jo Jo 
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thereby yielding SDE (4.5). 

(2) Next, suppose X satisfying SDE (4.5) is given. Since D is strictly in- 
creasing, X D i s _}_ = (X o D) s _ for any s > 0. Again, since any process is 
in synchronization with the continuous (^ r t)-time-change E, the 2nd changc-of- 
variable formula (3.1) applied to the integral form of SDE (4.5) yields 



(4.9) X t =x Q + p(E D{s _ h X D[s _)J)ds+ I a{E D ( s _),X D ( s _)J)dZ a 



Let Y X o D, then (4.9) immediately yields SDE (4.6), which completes the 



Remark 4.3. One may wonder whether the SDE 

dX t = p(Et,X t -)dt + ii(Et,X t -)dEt + <j(E t ,X t _)dZ Et with X Q = x 

can be reduced in the same manner as Theorem 4.2 (2). This is a question of 
whether the new driving process dt can be replaced by dDg t , which is possible 
only in very special cases; e.g., if D is continuous or p(E t , X t -) vanishes on 
every nonempty open interval D u ). 

For the remainder of this section, consideration mainly focuses on linear 
SDEs of the form 



(4.10) dX t = ( Pl (t, Et) + pa(*, E t )X t )dt + ( Ml (t, E t ) + p 2 (t, E t )X t )dE t 

+ (<T l (t,Et) + <T2(t,Et)Xt)dB Et with X = x a . 



Here B is a standard (^^-Brownian motion, E is a continuous (J-t)-time-change, 
and pj, pj, Oj (j = 1,2) are real- valued functions on R + xR + satisfying the 
following conditions: 



where L is a positive constant and Qt '■= J~E t - Note that a strong solution 
X to SDE (4.10) always has continuous paths due to the continuity of the 
driving processes. Conditions (4.3') and (4.4') respectively imply conditions 
(4.3) and (4.4); therefore, the uniqueness and existence of the strong solution 
X is guaranteed by Lemma 4.1. 

As demonstrated in the proof of Theorem 3.3, we have the handy calculus 
rules 




proof. 



□ 



(4.3') 
(4.4') 



\p 2 (t 7 u)\ + \p 2 (t 7 u)\ + \a 2 {t,u)\ < L for all f,«eM + , 
( Pj (t,E t )), {pj(t,E t )), (<Tj(t,Et)) G HQt) fOT j = 1,2 



(4.11) 



( 



[m, m] = [m, E] = [m, B o E] = [E, E] = [E, B o E) = 0, 
[B o E, B o E] = E, 
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where m denotes the identity map corresponding to Lebesgue measure. Remark 
4.3 implies that the simple substitution Y t ■= X Et fails to reduce even the most 
basic type of SDE (4.10) into a classical ltd SDE due to the presence of the dt 
term. This observation suggests that we establish a general form of solution to 
(4.10) via a direct approach rather than via such a simple substitution. It also 
calls into question the possibility of developing reduction schemes for converting 
SDEs of the form (4.10) into less complicated SDEs. Propositions 4.4 and 4.8 
together with Theorem 4.5 largely settle this issue. The linear SDE (4.10) is 
said to be homogeneous if p\ = \i\ = o\ = 0. 

Proposition 4.4. (Solution form for Homogeneous Linear SDEs) Let 

B be a standard {Tt)-Brownian motion. Let D and E be a pair satisfying 

[\D i — > E\\ or \[D < 1 E]]. Then the unique strong solution to the homogeneous 

linear SDE with initial condition 

(4.12) dX t = p 2 (t, E t )X t dt + n 2 (t, E t )X t dE t + a 2 (t, E t )X t dB Et , X Q = x Q 
is explicitly written as 

(4.13) X t =x cxpj^ p 2 (s,E s )ds + (jjn(8,E a )-^oi(a,E,))dE a 

+ o- 2 (s,E s )dB E ^, 

or equivalently as 

(4.14) X t = x expi.J p 2 (s,E s )ds + ([i 2 (D s -,s) -^a%(D s _,sj)ds 

+ * a 2 (D s _,s)dB s | . 

Proof. (4.14) follows from (4.13) together with the 2nd change-of- variable for- 
mula (3.1). Due to the uniqueness of the solution, it suffices to show that the 
process X given in (4.13) satisfies SDE (4.12). 

Let X be the process in (4.13) and write X t = x$ e At . A calculation similar 
to (3.10), via (4.11), yields [A, A] = o\{-,E) • E. By the Ito formula (1.1) with 
/(a) =x e a , 

(4.15) dX t = x e At dA t + \x e At d[A, A] t 

=X t {p 2 (t, E t )dt + (^(i, E t ) - i(T 2 2 (t, E t ))dE t + a 2 (t, E t )dB Et } 

+ \X t cj 2 2 {t,E t )dE t 
=p 2 (t, E t )X t dt + fi 2 (t, E t )X t dE t + a 2 (t, E t )X t dB Et . 

In addition, Xq = xq. Thus, X satisfies (4.12), completing the proof. □ 
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Theorem 4.5. (General Solution Form for Linear SDEs) Let B be a 

standard (Tt)-Brownian motion and D and E be a pair satisfying [[D i — > E\] 
or [[D < — i E}}. Then the unique strong solution to a general linear SDE (4.10) 
is explicitly written as 



(4.16) X t = *t 



x 



Pi{s,E s ) 
3, 



ds 



fj,i{s,E s ) - a 2 (s,E s )a 1 (s,E s ) 



dE, 



or equivalently as 



o~i(s,E s ) 



dB f 



(4.17) X t = $ t 



x 



/xi (D s - , s) - (72 (-Ds- , s)ai (D s - , s) 



ds 



DIs- 



3> 



_D(s- 



where $ is £/ie unique strong solution (4.13) to the homogeneous linear SDE 
(4.12) with xq replaced by 1. 3> is called the fundamental solution to the homo- 
geneous SDE (4.12). 

Proof. Since $ = 1 > 0, the explicit form (4.13) of $ shows that $ t > 
for all t > 0. Hence, the right hand side of (4.16) is meaningful. As in the 
proof of Proposition 4.4, it is sufficient to check that the process X in (4.16) 
satisfies SDE (4.10). For notational convenience, we suppress the dependence 
of the coefficients on E t and simply write pj(t) = pj(t,E t ), fJ-j(t) = fj,j(t,Et), 
a j (t)=a j (t,E t ) for j = 1,2. 

Let X be the process in (4.16) and write X t = &tZt- Since $ is the solution 
to SDE (4.12), the calculus rule (4.11) yields [$, Z] = (a 2 $ • (cti/$)) • E = 
(02CT1) • E. Hence, using the product formula (1.2), 

(4.18) dX t = $ t dZ t + Z t d<i> t + d[<S>,Z} t 

=pi{t)dt + (fjn(t) - a 2 (t)a 1 (t))dE t + a 1 (t)dB Et 

+ Zt(pa(t)*tdt + (i 2 (t)<PtdE t + a 2 {t)<$> t dB Et ) + o- 2 (t)ai(t)dE t , 

the right hand side of which yields that of SDE (4.10) upon replacing &tZ t by 
Xt. Moreover, Xq = xq, completing the proof. □ 

A multidimensional version of Theorem 4.5 can be obtained in a similar way 
by applying the Ito formula componentwise. 

Corollary 4.6. Let B be an n-dimensional standard [Tt)-Brownian motion 
starting at 0. Let (p 2 {t, E t )) , (fi 2 (t,E t )), (a 2 (t,E t )) (k = 1, . . . , ri) be dx d- 
matrix-valued processes, (pi(t,E t )), (ni(t,E t )), (a\(t,Et)) (k = l,...,n) be 
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d- dimensional processes. Let xq € M. d . Then the unique strong solution X to 
the SDE 

(4.19) dX t = (p 1 (t,E t )+p 2 (t,E t )X t )dt+(^ 1 {t,E t )+fX2(t,E t )X t )dE t 

n 

+ Y^{vx{t,Et)+(rl{t,E t )X t )dB k Et with X =x , 

k=l 

which is a d- dimensional process, is explicitly written as 

rt 



(4.20) 



•'•() 



+ f <i> s 1 p 1 {s,E s )ds 
Jo 

/ &J 1 (p^s, E s )-J2 E s )o- k i(s, E„))dE a 
J ° k=i 

/ K i y / ai(8,E B )dB% B 

JO u _, 



or equivalently as 
(4.21) X t = * t 



+ f $ s -Vi(s,S s )ds 
Jo 

^° fc=l 

Jo I — -I 



where $ = ($t) is t/ie fundamental solution to the homogeneous linear SDE 
corresponding to (4.19). Namely, $ is the unique d x d -matrix-valued process 
satisfying the homogeneous SDE 

n 

(4.22) d<Z> t =P2(t,Et)$tdt + p J2 (t,E t )$ t dE t + J24(t,Et)<f>tdBl ;t , 

k=l 

with initial condition <f>o = Id, where Id denotes the dxd -identity matrix. 

Proof. We first claim that for each path, <£>t is invertible for all t > 0. Otherwise, 
there would exist to > and A G R d \{0} such that $t A = 0. The d-dimensional 
process ($tA) satisfies the homogeneous linear SDE 

n 

(4.23) ^ t = / 9 2 (t,S t )* t * + M2(t,^t)*irf-Et+E (7 2(*,A)*td-B| t . 

fe=i 

The zero process is the unique solution to (4.23) for which ty to = £ IR d . 
Therefore, it follows that <&tA = for all t > 0, which contradicts $oA = A ^ 0. 
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Thus, <3? t is invertible for all t > 0, and the right hand side of SDE (4.20) is 
meaningful. 

As in the proof of Proposition 4.4, it suffices to show that X given in (4.20) 
satisfies SDE (4.19). Using the calculus rule 

[to, B k oE] = [E, B k oE\=Q and [B k oE,B e oE} = 5 kA E, 

where 5 k,i is the Kronecker delta, and applying the ltd formula componentwise, 
the proof is carried out in the same way as in Theorem 4.5. □ 

Remarks 4.7. (a) The advantage of rewriting solutions in the forms (4.14) 
and (4.17) is that they can be handled via conditioning on the random variable 
E t . This is especially useful in analyzing statistical data of a solution, such 
as its mean and variance. If Y is the solution to a classical Ito SDE with 
linear coefficients dY t = (b 1 (t) + b 2 (t)Y t )dt+ (n(t) + T 2 (t)Y t )dB t , then the first 
two moments of Y t are characterized as solutions to linear ordinary differential 
equations (ODEs), from which some information on statistics can be derived. 
(See [3, Thm. 4.5] for a general case, or [9, Problem 5.6.1] for a special case when 
t-2 = 0.) However, it is generally impossible to obtain such ODEs for the solution 
X to SDE (4.10), even when pj, pj, o~j are deterministic. For example, consider 
the SDE dX t = p 2 (t)XtdEt, a special case of (4.10). Taking expectations in 
the integral form, E[X t ] = xq +E[/ * p 2 {s)X s dE s \. The expectation and integral 
are not interchangeable due to the presence of the random integrator dE s . As a 
result, unlike the case of a classical Ito SDE, a general form of an ODE satisfied 
by E[A t ] cannot be obtained. This observation heightens the importance of 
expressions such as (4.14) and (4.17). Moreover, since these expressions are 
derived via the 2nd change-of- variable formula (3.1), there is no doubt that 
Formula (3.1) is an indispensible tool for dealing with SDEs of the form (4.1). 

(b) Recognizing how our new class of SDEs of the form (4.10) arise: View- 
point 1. If Et — t and pj, pj, o~j (j = 1,2) are all deterministic, then Proposition 
4.4 and Theorem 4.5 respectively reduce to well-known results for classical Ito 
SDEs with linear coefficients 

(4.24) dY t = (b 1 (t) + b 2 (t)Y t )dt+(a 1 (t)+a 2 (t)Y t )dB t , 

where bj(t) = pj(t) + pj(t) (j = 1,2). (See [3, Thm. 4.2].) This observation 
suggests that an SDE of the form (4.10) might be constructed via continuously 
altering the clock from t to Et in (4.24), but with the drift factor bj splitting 
into two components pj and pj, the former reflecting the effect of the original 
clock t and the latter of the new clock E t . Allocation of the weight of bj to pj 
and pj is due to consideration of how much the time-changed model is affected 
by the new clock. If the absolute value of pj is big (resp. small) in comparison 
to that of pj, then the model (4.10) contains a large (resp. small) effect of the 
original clock. SDE (4.5) with Z = B provides an example where there is no 
effect of the original clock. Note that pj and pj may take negative values as 
well. 
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Viewpoint 2. Again assume fj,j, o~j (j = 1,2) are all deterministic. Adopt a 
classical ltd SDE 

(4.25) dZ t = ( Pl (t) + p 2 {t)Z t )dt + ( Tl (t) + T 2 (t)Z t )dB t 

as the starting form of SDE (4.10). This interpretation is valid when path 
properties or statistical data of the solution to a simple SDE of the form (4.25) 
fail to match the real data (especially in terms of the volatility coefficients Tj), 
but clearly possesses a drift similar to (pi(t) + p2(t)Z t )dt. In this situation, 
one prefers to 'break' the dB t term via changing the clock from t to E t so that 
the model has more flexibility in describing the volatility. As a result, dE t and 
dB Et terms are obtained as in (4.10), without changing the drift coefficients pj. 
Note that the arguments from both viewpoints apply to a general class of SDEs 
of the form (4.1) as well. 

(c) The general form of solutions obtained in Proposition 4.4, Theorem 4.5 
and Corollary 4.6 are all valid even when SDE (4.10) has general process coef- 
ficients. More precisely, if the coefficients pj, fj,j, dj (j = 1,2) are processes in 
L(<7t) with Qt '■= J-E t such that the absolute values of p 2 , /i2, o~ 2 arc dominated 
by some random variable L, then it can be shown by reformulating Theorem 7 
of [13, Chap. V] that SDE (4.10), with the coefficients evaluated at (t,co) rather 
than (t, E t (uj)) , has a unique strong solution; moreover, the explicit form of the 
solution has exactly the same expression as in the previous results. 

Just as there is a reduction method for classical Ito SDEs with nonlinear 
coefficients 

(4.26) dY t = b(t, Y t )dt + T {t)Y t dB t with Y = x , 

Proposition 4.8 provides an analogous technique for approaching a certain type 
of nonlinear SDE including terms driven by a time-changed Brownian motion. 
The 'integrating factor' 

U t :=exp{^J T 2 (s)ds- j T(s)dB s } 

reduces (4.26) to a path-by-path ODE d(U t Y t ) = U t ■ b(t, Y t )dt, with U Y = x , 
computation of which almost traces the proof of Proposition 4.8. Applications 
of this reduction scheme arc provided in Examples 5.5 and 5.6. 

Proposition 4.8. (Reduction Method) Let B be a standard (Tt)- Brownian 
motion. Let E be a continuous (J-t) -time- change. Then the 'integrating factor' 
U defined by 

(4.27) C/ t :=cxpj^ (^o-l{s,E s ) - p 2 (s,E s ))dE a - a 2 {s, E s )dB Es j 
reduces the nonlinear SDE 

(4.28) dX t = p(t, Et,X t )dt + p 2 (t, E t )X t dE t + a 2 (t, E t )X t dB Et , A = x , 
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to a path-by-path ODE 

(4.29) — - = U t - P (t,E u U^W t ), Wo = x , 



where W t := U t X t so that X t = U^Wf 

Proof. For notational convenience, we suppress the dependence on E t and sim- 
ply write p(t,X t ) = p(t,E t ,X t ), p 2 {t) = fJ>2(t,E t ) and a 2 (t) = a 2 (t,E t ). Write 
U t = e At so that 

A t = J (^of (s) - p 2 (s))dE s - J a 2 {s)dB Es . 

Then the Ito formula (1.1) with f(a) = e a together with the calculus rules in 
(4.11) yield 

dU t = U t dAt + \U t d[A,A} t = U t {(<4(t) - p 2 {t))dE t - a 2 (t)dB Et }. 

Hence, by the product formula (1.2), 

(4.30) d(U t X t ) = U t dX t + X t dU t + d[U, X] t 

= U t { P (t, X t )dt + p 2 (t)X t dE t + <j 2 (t)X t dB Et } 

+ X t U t {(4(t) - p 2 {t))dE t - a 2 {t)dB Et ) - a 2 2 (t)X t U t dE t 
= U t ■ p[t, X t )dt. 

By setting W t := U t X t , (4.30) immediately yields (4.29). □ 

5 Examples 

The examples below are drawn from the classical Ito SDEs; however, the driving 
processes involve a continuous time-change E and the time-changed Brownian 
motion B o E. Assume that all coefficients of SDEs appearing in this section 
satisfy the conditions (4.3) and (4.4). 

Example 5.1. The most basic linear SDE is the homogeneous one with constant 
coefficients, which is an analogue of the so-called Black-Scholcs SDE. Consider 

(5.1) dX t = pX t dt + pX t dE t + oX t dB Et with X = x , 

where p, p. a are real constants and x > 0, a > 0. The case where E t = t 
corresponds to the Black-Scholes model dY t = (p + p)Y t dt + aY t dB t , which 
describes the price Y t of a risky asset — say a stock price — with an expected 
rate of return b '■= p + p. For example, let us assume b > 0. Then the current 
stock price being positive (resp. negative) means that the drift coefficient bY t 
is also positive (resp. negative); hence, the stock price Y t is expected to drift 
in the positive (resp. negative) direction. The other coefficient a exhibits the 
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variability of the stock model which is also proportional to the current level of 
the price. 

In the classical Black-Scholcs SDE dY t = bY t dt + aY t dB t with Yq = xo, its 
solution 



Y t = x exp{(&- ^<7 2 y + aB t } 



has the following asymptotic behavior: 
(Y.l) If b > cr 2 /2, then limt^oo Y t = oo. 
(Y.2) If b < a 2 /2, then lim^oo Y t = 0+. 

(Y.3) If b — er 2 /2, then Y t asymptotically fluctuates between arbitrarily large 
and arbitrarily small positive values infinitely often. 

This follows by rewriting the solution as Y t = xq cxpji \(b — ct 2 /2) + a ■ B t /t\ } 
and using the law of the iterated logarithm for paths of Brownian motion 



(5.2) limsup , = = 

t->oo V2* log log* 



1 and lim inf ■ 



Bi 



t^OO y/2t lOg lOg t 



Analysis of the asymptotic behavior of the solution to SDE (5.1) is accom- 
plished with the help of the explicit solution form obtained from (4.13), 



(5.3) 



X t = x a cxpjpt + Ui- ^er 2 ) E t + aB Et | . 



First, if p = 0, i.e., if there is no effect of the original clock upon the solution 
of SDE (5.1), then, by Theorem 4.2 (2), (X Dt ) satisfies the classical Ito SDE 
dX[) t = fiX£) t dt+aXi) t dB t . Since lim^oo D t ~ oo, X has the same asymptotic 
behavior as the above-mentioned Y with b replaced by \x: 

(X.a.l) If p = and p, > cr 2 /2, then lim^oo X t = oo. 

(X.a.2) If p = and p < cr 2 /2, then lim^oo X t = 0+. 

(X.a.3) If p — and p = <r 2 /2, then X t asymptotically fluctuates between 
arbitrarily large and arbitrarily small positive values infinitely often. 

Next, suppose p ^ 0. Assume lim^oo E t — oo and lim^oo E t /t = 0; i.e., 
E t is asymptotically slower than t. By rewriting (5.3) as 



X t = .T exp<^ t 



1 



V- 7; a -r 



E, 



B 



E t 



E, 



and using (5.2) again, we easily observe 



(X.b.l) If p > and E t is asymptotically slower than t, then lim^oo X t 
(X.b.2) If p < and E t is asymptotically slower than t, then limt-.oo X t 



oo. 
0+. 
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These cases match with our intuition: if the original clock t asymptotically ticks 
more frequently than the new clock E t , then the p describing the effect of the 
original clock completely determines the future behavior of the solution X, no 
matter what values p and a take. 

On the other hand, if E t grows faster than t, i.e., if lim^oo E t /t — oo, then 
the situation becomes much more complicated. Rewrite (5.3) as 



X t = x Q exp< E t 



t 



1 r 



.Be, 
E t 



By noting (5.2) again, we observe 

(X.c.l) If p ^ 0, p > er 2 /2 and E t grows faster than t, then lim^oo X t = oo. 

(X.c.2) If p ^ 0, p < er 2 /2 and E t grows faster than t, then lim^oo X t = 0+. 

(X.c.3) If p ^ 0, p = cr 2 /2 and Et grows faster than t, then the fluctuation of 
Xt varies depending on the coefficients of the SDE and also the speed 
at which E t grows. 

The first two cases show that if p ^ a 2 / '2 and E t grows faster than t, then 
the asymptotic behavior of X, regardless of the value of p(^ 0), coincides with 
(X.a.l) and (X.a.2). This is due to the fact that the effect of the faster clock E t 
is strongly reflected on p to the extent that p is ignored. 

In the special situation (X.c.3), if limt_>oo \/2Et log log Et/t = oo so that 
E t grows extremely fast, then X t asymptotically takes arbitrary values on the 
positive real line infinitely many times. This is immediate upon writing 



X t = x cxp<! t 



D 



^2E t \og\ogE t 



^2E t \og\og E t 
t 



On the other hand, if, e.g., lim^oo \/2E t log log E t /t = 0, then X t asymptoti- 
cally goes off to oo if p > and decreases to if p < 0. 

These observations establish that as the time-change E accelerates the speed 
at which time passes, dependence of the behavior of the solution X upon p and 
p respectively becomes lighter and heavier. □ 

Example 5.2. Assume B is independent of D, or equivalently, of E. The 
homogeneous linear SDE 



(5.4) 



dX t = p{t)X t dt + p(E t )X t dE t + o{E t )X t dB Et with X = x , 



where xq > 0, has a unique strong solution X expressed as (4.14). 

The value of the mean function E[X t ] can be investigated by conditioning 



24 



on E t and using the independence of B and E: 



(5.5) E[X t ] = x Q cxpj / p(s)ds| 



x E 



exp 



xq exp 



cxp{ ^ ( M (s) - -a 2 (s))ds + jf <r(s)dB a } 
p(s)ds| 

f (MW-^ 2 (s))ds + jT ct(s)cLB s }] pt(d«) 
/i(s)dsj -E[M„] Pt(efo), 



o 

exp 
/ 



where pt denotes the law of the random variable E t and M is a continuous 
(7 c t)-local martingale given by 



(5.6) 



M v := exp< — 



HP 



(s)ds 



i 

a(s)dB s y 



Actually M is a martingale since er satisfies the Novikov condition; i.e., 
E[exp{i/ "<T 2 (s)ds}] < oo for all w > 0. (See [9, Prop. 3.5.12].) Hence, 
E[M V ] = 1 for all v > 0. Thus, (5.5) yields 

(5.7) E[X t ] = xo exp|y p(s)ds| • y exp p,(s)ds^ pt(dv). 

If E't = t, then p< = <5t, the Dirac measure with mass at t. Hence, (5.7) yields 
E[Xj] = a;oexp{J (p(s) + p(s)^)ds}, which, of course, coincides with the mean 
function E[Y t ] of the solution Y to the classical Ito SDE dY t = (p{t)+p(t))Y t dt+ 
a(t)Y t dB t with Y = xq. The result (5.7) shows that the behaviors of p and p 
together govern the range of fluctuation of the mean function E[Xt]. Moreover, 
even when the coefficient of dEt in SDE (5.4) is replaced by a more general 
p(t, Et)Xt, some form of estimate on M[Xt] can still be obtained. For instance, 
if L p(D s _,s)ds > for all v > 0, then E[X t ] > x cxp{f Q p(s)ds\. 

The variance function VLY t ] of the solution X is computed similarly: 



V[X t ] =x 2 exp{2 / p(s)ds} 



exp 



|2 J p(s)ds + J o 2 (s)ds^ p t {dv) 
H(s)dsj p t (dv) J 



exp 



Unlike the explicit form of the mean function in (5.7), V[X*] involves the infor- 
mation a concerning the weight of the dBE t term in SDE (5.4). 
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As a special case of SDE (5.4), assume fj,(u) = —A for some A > 0. Then 

(5.7) is expressed in terms of the Laplace transform of the law of E t : 

(5.8) E[X t ] = x Q exp{^ p{s)ds} ■ J e~ Xv Pt (dv). 

Moreover, if E is the first hitting time process of an (jT t )-stable subordinator of 
index (3 £ (0, 1) which is independent of B, then the Laplace transform in (5.8) 
is associated with the Mittag-Lcffler function due to [2, Thm. 4.3]: 

(5.9) E[X t ] = x cxp{^ P (s)ds} • Ep(-\t% 

where E^^) := X^S^Lo z n jT{fin + 1) with L(-) being the Gamma function. □ 
Example 5.3. Consider the inhomogeneous linear SDE 

(5.10) dX t = (J^ - Y L -X t )dt + - ^tfdEk + dB*, 

t £ [0, 1) with X Q = a, 

where a, 6, c, 7, rj £ K and E t increases to 1 as t increases to 1. 

The fundamental solution to the homogeneous linear SDE corresponding to 
(5.10) is $ t = (1 - ty<{\ - Etf. Hence, (4.16) yields 

,,n, ^ (1 _ tr(1 _ w „ + ^ ( l^ ( ^f)» d8 
+ / t^:(— ) i—J 

If Et = t, then the solution (5.11) reduces to 

r* 1 

(5.12) X t = (l-t)a-t(b + c) + (l-t) dB„ 

io l-« 

a Brownian bridge from a to (6 + c). Moreover, the class of SDEs of the form 
(5.10) contains a 'time-changed Brownian bridge' from a to c. In fact, if b = 0, 
7 = and rj = 1, then X satisfies the SDE 

(5.13) dX t = f—^-r - —^—X t )dE t + dB Et , t £ [0, 1) with X = a, 

\ 1 — E t 1 — E t I 

which is, by Theorem 4.2, associated with the classical Brownian bridge SDE 

(5.14) dY t = (j^— - —Y^jdt + dBt ' 1 e t°> l ) with X ° = a - 

via the relation X = YoE. Thus, in this particular case, X is a process obtained 
by time-changing the Brownian bridge Y. □ 



26 



Viewpoint 1 of Remark 4.7 (b) states that it is possible to recognize that 
the two components pj and fij of SDE (4.10) are produced by splitting the drift 
factor of some classical Ito SDE. Examples 5.1, 5.2 and 5.3 are all discussed 
from this viewpoint. However, as mentioned in Viewpoint 2 of the remark, it is 
also possible to attribute the presence of \ij to the dB t term in a classical Ito 
SDE. Example 5.4 illustrates this viewpoint. 

Example 5.4. This example investigates statistical data obtained from the 
solution to the inhomogencous linear SDE 

(5.15) dX t = —aX t dt + [idEt + adBE t with Xq = Xq, 

where a, a > 0, \i £ R, and xq ^ 0. SDE (5.15) with E t = t and /j, = is 
called the Langcvin equation or the Ornstein-Uhlcnbcck model, and its solution 
is referred to as the Ornstein-Uhlenbeck process. The coefficient — aX t of the 
dt term is negative (resp. positive) when X t is positive (resp. negative), which 
implies X t is drawn back to zero once it drifts away. Since the coefficient \x 
describing the drift based on the new clock E t is not proportional to the current 
position Xt, if, e.g., E t represents the business time at the calendar time t, then 
Xt, regardless of its value, is always affected by the evolution of the business 
time. In other words, the model has a certain factor of weight [i which pushes the 
position either up or down during business hours, and its effect on the position 
becomes larger (resp. smaller) when the business time grows faster (resp. slower). 
Moreover, the dispersion coefficient a does not depend on the position cither. 
Therefore, unless the time-change E either stays flat or accelerates or decelerates 
drastically on an interval, X t fluctuates on this interval at a certain rate with 
mild error even when it approaches close to zero. In finance, the Ornstcin- 
Uhlcnbeck-type model (5.15), which incorporates a possible time-change, could 
be used to describe the deviation of an interest rate around a central bank's 
target rate. 

Assume both of the following technical conditions are satisfied: 

(a) each random variable E t is bounded; i.e., f(E t < Ct) = 1 for some finite 
positive constant Ct; 

(b) E[/ * e 2aD '-ds] < oo for all t > 0. 

The monotonicity of D implies that the condition (b) is equivalent to: 

(b') E[e 2aDt -] < oo for all t > 0. 

Let us analyze the mean ~E[X t ] of the solution X to SDE (5.15). By (4.16) 
and (4.17), X can be represented in two ways: 

(5.16) X t = e- at {x +iJ,J e as dE s +aJ^ e as dB Es } 



e -at. 
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By assumption (b), the process TV defined by N t ■= J Q e aD "~dB s is an {T t )- 
martingale. Since each E t is a bounded (^ r t )-stopping time due to (a), Doob's 
optional sampling theorem yields E[NE t ] = E[N ] = 0. (Sec [9, Problem 1.3.23 
(i)].) Hence, taking expectations in (5.16), 



(5.17) 



E[X t ] =e' at { x + Li 



e as dE s 



z aD °-ds 



Consequently, the asymptotic behavior of the mean function E[X t ] completely 
depends on the distributions of the processes E and D. In the special case 
where E t (uj) = R(uj) ■ t for some positive random variable R, E[X t ] = x n e~ at + 
(fj,E[R]/a)(l — e~ at ), which approaches iM[R]/a as t — > oo. Therefore, if the 
force attracting X t to zero is sufficiently strong compared to the factor producing 
the effect of the evolution of the time (i.e., if a is much larger than the absolute 
value of /i and E[i?]), then the expected value of the position tends to a level 
close to zero as t — > oo. On the other hand, the bigger the weight /i or the 
expected rate E[i?] of acceleration of the new clock, the greater the asymptotic 
value of the expected position. 

Another way to observe the fluctuation of E[JQ] is to directly analyze the 
integral form of the SDE (5.15). Taking the expectation, 



(5.18) 



E[X t ] = -a [ E[X s ]ds + (jE[Et]+aE[B Et ]. 
Jo 



The last term vanishes again due to the assumption (a) and Doob's optional 
sampling theorem. Hence, we obtain a differential equation 

(5.19) ^E[X t ] = -aE[X t ] + fi^-E[E t ] with E[X ] = x , E[E ] = 0. 
dt dt 

Although this is not the explicit form of E[X t ] obtained in (5.17), it still provides 
information on the relationship between the time evolutions of E[Xt] and E[_Ej]. 

The term E[2?E t ] in (5.18) vanishes even when the assumption (a) is replaced 
by one of the following: 

(c) EyW t ] < oo for all t > 0; 

(d) B is independent of E. 

If condition (c) holds, which is weaker than (a), then the 'Wald identity' E[Be 4 ] 
= holds for each t > 0. (See [9, Problem 3.2.12, Exercise 3.3.35].) On the 
other hand, (d) encourages conditioning on the random variable E t to obtain 
E[B Et ] = 0. 

Suppose E is the first hitting time process of an (J r t )-stable subordinator of 
index (3 € (0, 1) which is independent of B, so condition (d) holds by assumption. 
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There is a positive constant c(f3) such that K[E t ] = c(/3) t 13 for all t > 0, due to 
[10, Cor. 3.1]. Hence, (c) also holds. Moreover, using this moment result, (5.19) 
is reexprcsscd as 



(5.20) ±E[X t ] 



aE[X t }+ fiPciP)^- 1 with E[X ]=x . 



The solution of the first order linear ODE (5.20) is given by 



(5.21) 




Jo 



where g a ,t(r) '■= e a ^~ r \ and T(-) and J 13 respectively denote the Gamma func- 
tion and the fractional integral of order [3. (For definition of fractional integrals, 
see [4].) 

An interesting conjecture can be made by comparing SDE (3.16) in Example 
3.6 and SDE (5.15), both for the particular E discussed in the above paragraph 
which is assumed independent of B. First, SDE (5.15) is particularly different 
from SDE (3.16) due to the presence of the dt term. Second, Theorem 4.1 in 
Hahn, Kobayashi and Umarov [5] shows that the transition probability density 
of the solution to SDE (3.16) satisfies PDE (3.17), and the proof is carried out 
by taking the expectation in the time-changed Ito formula (3.5). Consequently, 

(5.21) suggests that if SDE (3.16) is replaced by an SDE having a term p(X t )dt, 
then the corresponding PDE may involve a fractional integral term. □ 

The following two examples clarify how to apply the reduction method ob- 
tained in Proposition 4.8. 

Example 5.5. Solution (5.3) to the homogeneous linear SDE (5.1) discussed in 
Example 5.1 can also be obtained by using the technique provided in Proposition 
4.8. In this case, the integrating factor is Ut = exp{(<7 2 /2 — p,)E t — aBE t } and 
(4.29) becomes the path-by-path ODE dW t = pW t dt with W = x , which has 
the solution W t — Xoe pt . Hence, the relation X t = U^Wt immediately yields 
the desired solution form (5.3). More generally, the same reduction scheme 
proves Proposition 4.4. □ 

Example 5.6. As another application of the reduction method introduced in 
Proposition 4.8, consider a generalized population growth model 

(5.22) dX t = qX t (K - X t )dt + pX t dE t + <rX t dB Et with X = x 

where q, K, xq > and /x, a G K. This model describes the growth of a 
population of size X t in some environment, q and K represent the quality and 
the carrying capacity of the environment, respectively. If the quality of life is 
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good and the current population is less than the carrying capacity, i.e., if q is 
large and < X t < K, then the population will grow, i.e., the drift coefficient 
qX t (K—X t ) is positive. On the other hand, a population exceeding the capacity 
of the environment is expected to decrease even when the quality is good, i.e. 
if X t > K, then the drift qX t (K — X t ) is negative, regardless of the value of 
<Z(>0). 

Note that SDE (5.22) possesses a distinct form of coefficients in dt and dE t 
terms, unlike Examples 5.1, 5.2 and 5.3. Hence, this model is constructed based 
on Viewpoint 2 of Remark 4.7 (b). The presence of the term /j,X t dE t implies that 
a certain factor originating in the new clock affects the growth of the population, 
and the effect is proportional to the current position X t . o describes the noise 
of the system as in the classical population growth model (i.e., SDE (5.22) with 
E t =t and — 0). 

Theorem 4.5 cannot be applied to the nonlinear SDE (5.22). Instead, Propo- 
sition 4.8 with W t — U t X t where U t = exp{(er 2 /2 — fi)E t — aB Et }, yields the 
path-by-path ODE 

(5.23) ^ = qWAK - U^Wt) with W a = x . 

dt 

Consider a Bernoulli-type ODE 

(5.24) y'(t) = f(t)y 2 (t) + ky(t) with y(0) = x , 

where A: is a real constant and the symbol ' denotes the derivative with respect 
to t. By the substitution z(t) = y _1 (i), the ODE (5.24) reduces to z'(t) + 
kz(t) = —f(t) with z(0) = Xq . Multiplication of both sides by e kt leads to 
{e ht z(t)Y = —e kt f(t), whose solution is 

rt kt 

(5.25) e kt z(t)^x^ = - e ks f( S )ds, or y(t) = — — . 

Jo Xq J e j (s)ds 

The substitutions, y{t) = W t , f(t) = -qll^ 1 , k = qK in (5.25), yield 

Uf 1 ■ exp{qKt} 



(5.26) X t = U t ~ L Wt 



x o 1 + Jo G^p{qKs} ■ qU s 1 ds 

exp{qKt + ((J,- \a 2 )E t + aB Et } 
x o 1 + 1 Jo expjgif s + (/i — \a 2 )E s + oB Es }ds 

□ 



Appendix — Construction of Stochastic Integrals Driven 
by a Semimartingale 

The aim of this appendix is to make explicit the class L(Z,Tt) of Z-integrable 
predictable processes treated in this paper. For details regarding the construc- 
tion of stochastic integrals, consult [13, II-IV]. 
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Throughout, a filtration (Ft) satisfying the usual conditions is fixed. Write 
D = V)(T t ) (cadlag adapted processes), L = h(T t ) (caglad adapted processes), 
and V = V{Tt) (predictable processes). Let bL and \>V denote bounded pro- 
cesses in the specified class. Let § be a subset of L consisting of all processes of 
the form H t = H Q l^(t) + Y?i=i HiI(T it T i+1 ](t), where n is a positive integer, 
{Ti}™^ 1 is an increasing sequence of finite stopping times with T\ = 0, and each 
Hi is an Tt± -measurable random variable. 

First, endow KD, L and § with the topology induced by "H m — ► H if and 
only if for each t > 0, sup 0<s<t \H™ — H s \ — > in probability aswn oo." 
Then § is a dense subspace of L, and D becomes a complete metric space 
with a compatible metric d(Y, Z) := ][]^ =1 (l/2 ,l )E[min(l, sup < s < t \Y S - Z s \)] . 
Given a semimartingale Z starting at 0, the stochastic integral of H G § of the 
above form is defined to be H • Z := J Z (H) := Hi(Z T '+ l - Z T *) where 

Z^ := Z m [ n ( t T y The continuous linear operator Jz ■ § — > D uniquely extends 
to an operator defined on L. For the moment, denote Jz(H) as [Dl-]iJ • Z 
for H G L. Note that the quadratic variation of Z is defined by (2.1) via this 
integral operator. 

The next step is to introduce the^ space H 2 ofjBemimartingales starting at 
with a unique decomposition Z = M + A where M is a local martingale and A 
is a predictable process of finite variation such that 



\Z\ 



\[M,M]V?\ 



L- 



\dA s \ \\ L 2 < oo. 



The real vector space H 2 with the norm || • ||^2 forms a Banach space. To extend 
a class of integrands, first fix an integrator Z = M + A G TL 2 and introduce a 
metric d^ on hV by 



(H s -K s ) 2 d[M,M} s j ^+ J \H s -K s \\dA 



where |gL4 s | denotes the integral with respect to the total variation measure. 
The integrals appearing in this definition are understood path-by-path in the 
Lebesgue-Stieltjes sense, and it follows that d^(H, K) = \\H • Z — K • Z\\^2. 
Under this metric, bL is dense in hV. For H G b"P, it is easy to see that a 
unique 7i 2 -limit of the sequence {\Dl-]H n • Z} exists where {H n } is an approx- 
imating sequence in bL for H. Moreover, the limit is determined independently 
of the choice of the approximating sequence. Hence, the stochastic integral 
[D2]-H •Z := H 2 - Iim n ^ 00 [Dl-].ff' n • Z is well-defined. 

The third step requires another class of integrands, denoted L^i^Z^Tt), 
which consists of predictable processes with 



H 



'■d[M,M] s } 



1/2 



+ 


/>oo 

/ \H s \\dA s \ 


< oo 


L 2 


Jo 


L 2 



where Z = M + A G TL 2 . Associate to H G L u -2 (Z, JF t ), the truncation processes 
{H k } in hV, given by H h := H\s\m<k}- Again, via the same reasoning as 



31 



above, the stochastic integral [D3-]H • Z is defined to be the unique H 2 -limit of 
the sequence {[D2-]H k • Z}. That is, [D3]-# • Z := H 2 - lim„^ 00 [D2-]i?" • Z. 

Finally, given a general scmimartingale Z starting at 0, a predictable process 
H is said to be Z -integrable, denoted H £ L(Z,Tt), if there exists a sequence 
{<r n } of stopping times increasing to oo such that Z" := Z°""~ £ 7i 2 and 
H £ L H *(Z n ,F t ) for each n, where Z t a ~ := Z t I [0j(T) (i) + Z CT _ I [<T)00 )(f). With 
this sequence {cr™}, the stochastic integral of if driven by Z is defined to be 
H • Z .= [D3-]iJ • Z" on [0, er n ). This definition is consistent and independent 
of the choice of the localizing sequence {c™}. 

One important special case is when Z = M is a continuous (jF t )-local martin- 
gale. In this case, H £ L(M, T t ) if and only if H £ P(JF { ) and P( /„' H 2 d\M, M] s 
< oo) = 1 for all t > 0. Moreover, the stochastic integral H»M is also a contin- 
uous (J-t)-local martingale. In particular, if Z = B is a standard (J r t )-Brownian 
motion and E is a continuous (^ r t)-time-change, then it is easily shown that 
(Be,.) is a continuous ((Jt)-local martingale, where Qt '■= J~E t - Thus, for any 
K £ L(B o E,Q t ), the stochastic integral K • (B o E) is also a continuous 
(£7t)-local martingale. 
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